Long-term experimental evolution in Escherichia coli. X. Quantifying the fundamental and realized niche

BACKGROUND: Twelve populations of the bacterium, Escherichia coli, adapted to a simple, glucose-limited, laboratory environment over 10,000 generations. As a consequence, these populations tended to lose functionality on alternative resources. I examined whether these populations in turn became inferior competitors in four alternative environments. These experiments are among the first to quantify and compare dimensions of the fundamental and realized niches. RESULTS: Three clones were isolated from each of the twelve populations after 10,000 generations of evolution. Direct competition between these clones and the ancestor in the selective environment revealed average fitness improvements of ~50%. When grown in the wells of Biolog plates, however, evolved clones grew 25% worse on average than the ancestor on a variety of different carbon sources. Next, I competed each evolved population versus the ancestor in four foreign environments (10-fold higher and lower glucose concentration, added bile salts, and dilute LB media). Surprisingly, nearly all populations were more fit than the ancestor in each foreign environment, though the margin of improvement was least in the most different environment. Most populations also evolved increased sensitivity to novobiocin. CONCLUSIONS: Reduced functionality on numerous carbon sources suggested that the fundamental niche of twelve E. coli populations had narrowed after adapting to a specific laboratory environment. However, in spite of these results, the same populations were competitively superior in four novel environments. These findings suggest that adaptation to certain dimensions of the environment may compensate for other functional losses and apparently enhance the realized niche

Evolutionary rates and gene dispensability associate with replication timing in the Archaeon Sulfolobus islandicus

In bacterial chromosomes, the position of a gene relative to the single origin of replication generally reflects its replication timing, how often it is expressed, and consequently, its rate of evolution. However, because some archaeal genomes contain multiple origins of replication, bias in gene dosage caused by delayed replication should be minimized and hence the substitution rate of genes should associate less with chromosome position. To test this hypothesis, six archaeal genomes from the genus Sulfolobus containing three origins of replication were selected, conserved orthologs were identified, and the evolutionary rates (dN and dS) of these orthologs were quantified. Ortholog families were grouped by their consensus position and designated by their proximity to one of the three origins (O1, O2, O3). Conserved orthologs were concentrated near the origins and most variation in genome content occurred distant from the origins. Linear regressions of both synonymous and nonsynonymous substitution rates on distance from replication origins were significantly positive, the rates being greatest in the region furthest from any of the origins and slowest among genes near the origins. Genes near O1 also evolved faster than those near O2 and O3, which suggest that this origin may fire later in the cell cycle. Increased evolutionary rates and gene dispensability are strongly associated with reduced gene expression caused in part by reduced gene dosage during the cell cycle. Therefore, in this genus of Archaea as well as in many Bacteria, evolutionary rates and variation in genome content associate with replication timing

Why genes evolve faster on secondary chromosomes in bacteria

In bacterial genomes composed of more than one chromosome, one replicon is typically larger, harbors more essential genes than the others, and is considered primary. The greater variability of secondary chromosomes among related taxa has led to the theory that they serve as an accessory genome for specific niches or conditions. By this rationale, purifying selection should be weaker on genes on secondary chromosomes because of their reduced necessity or usage. To test this hypothesis we selected bacterial genomes composed of multiple chromosomes from two genera, Burkholderia and Vibrio, and quantified the evolutionary rates (dN and dS) of all orthologs within each genus. Both evolutionary rate parameters were faster among orthologs found on secondary chromosomes than those on the primary chromosome. Further, in every bacterial genome with multiple chromosomes that we studied, genes on secondary chromosomes exhibited significantly weaker codon usage bias than those on primary chromosomes. Faster evolution and reduced codon bias could in turn result from global effects of chromosome position, as genes on secondary chromosomes experience reduced dosage and expression due to their delayed replication, or selection on specific gene attributes. These alternatives were evaluated using orthologs common to genomes with multiple chromosomes and genomes with single chromosomes. Analysis of these ortholog sets suggested that inherently fast-evolving genes tend to be sorted to secondary chromosomes when they arise; however, prolonged evolution on a secondary chromosome further accelerated substitution rates. In summary, secondary chromosomes in bacteria are evolutionary test beds where genes are weakly preserved and evolve more rapidly, likely because they are used less frequently

Comparative genomics of Burkholderia multivorans, a ubiquitous pathogen with a highly conserved genomic structure

The natural environment serves as a reservoir of opportunistic pathogens. A well-established method for studying the epidemiology of such opportunists is multilocus sequence typing, which in many cases has defined strains predisposed to causing infection. Burkholderia multivorans is an important pathogen in people with cystic fibrosis (CF) and its epidemiology suggests that strains are acquired from non-human sources such as the natural environment. This raises the central question of whether the isolation source (CF or environment) or the multilocus sequence type (ST) of B. multivorans better predicts their genomic content and functionality. We identified four pairs of B. multivorans isolates, representing distinct STs and consisting of one CF and one environmental isolate each. All genomes were sequenced using the PacBio SMRT sequencing technology, which resulted in eight high-quality B. multivorans genome assemblies. The present study demonstrated that the genomic structure of the examined B. multivorans STs is highly conserved and that the B. multivorans genomic lineages are defined by their ST. Orthologous protein families were not uniformly distributed among chromosomes, with core orthologs being enriched on the primary chromosome and ST-specific orthologs being enriched on the second and third chromosome. The ST-specific orthologs were enriched in genes involved in defense mechanisms and secondary metabolism, corroborating the strain-specificity of these virulence characteristics. Finally, the same B. multivorans genomic lineages occur in both CF and environmental samples and on different continents, demonstrating their ubiquity and evolutionary persistence

Bacterial community profiles and Vibrio parahaemolyticus abundance in individual oysters and their association with estuarine ecology

Oysters naturally harbor the human gastric pathogen Vibrio parahaemolyticus, but the nature of this association is unknown. Because microbial interactions could influence the accumulation of V. parahaemolyticus in oysters, we investigated the composition of the microbiome in water and oysters at two ecologically unique sites in the Great Bay Estuary, New Hampshire using 16s rRNA profiling. We then evaluated correlations between bacteria inhabiting the oyster with V. parahaemolyticus abundance quantified using a most probable number (MPN) analysis. Even though oysters filter-feed, their microbiomes were not a direct snapshot of the bacterial community in overlaying water, suggesting they selectively accumulate some bacterial phyla. The microbiome of individual oysters harvested more centrally in the bay were relatively more similar to each other and had fewer unique phylotypes, but overall more taxonomic and metabolic diversity, than the microbiomes from tributary-harvested oysters that were individually more variable with lower taxonomic and metabolic diversity. Oysters harvested from the same location varied in V. parahaemolyticus abundance, with the highest abundance oysters collected from one location. This study, which to our knowledge is the first of its kind to evaluate associations of V. parahaemolyticus abundance with members of individual oyster microbiomes, implies that sufficient sampling and depth of sequencing may reveal microbiome members that could impact V. parahaemolyticus abundance

Genetic characterization of clinical and environmental Vibrio parahaemolyticus from the Northeast USA reveals emerging resident and non-indigenous pathogen lineages

Gastric infections caused by the environmentally transmitted pathogen, Vibrio parahaemolyticus, have increased over the last two decades, including in many parts of the United States (US). However, until recently, infections linked to shellfish from the cool northeastern US waters were rare. Cases have risen in the Northeast, consistent with changes in local V. parahaemolyticus populations toward greater abundance or a shift in constituent pathogens. We examined 94 clinical isolates from a period of increasing disease in the region and compared them to 200 environmental counterparts to identify resident and non-indigenous lineages and to gain insight into the emergence of pathogenic types. Genotyping and multi-locus sequence analysis (MLSA) of clinical isolates collected from 2010 to 2013 in Massachusetts, New Hampshire, and Maine revealed their polyphyletic nature. Although 80% of the clinical isolates harbored the trh hemolysin either alone or with tdh, and were urease positive, 14% harbored neither hemolysin exposing a limitation for these traits in pathogen detection. Resident sequence type (ST) 631 strains caused seven infections, and show a relatively recent history of recombination with other clinical and environmental lineages present in the region. ST34 and ST674 strains were each linked to a single infection and these strain types were also identified from the environment as isolates harboring hemolysin genes. Forty-two ST36 isolates were identified from the clinical collection, consistent with reports that this strain type caused a rise in regional infections starting in 2012. Whole-genome phylogenies that included three ST36 outbreak isolates traced to at least two local sources demonstrated that the US Atlantic coastal population of this strain type was indeed derived from the Pacific population. This study lays the foundation for understanding dynamics within natural populations associated with emergence and invasion of pathogenic strain types in the region

Phylogenomic study of Burkholderia glathei-like organisms, proposal of 13 novel Burkholderia species and emended descriptions of Burkholderia sordidicola, Burkholderia zhejiangensis, and Burkholderia grimmiae

Partial gyrB gene sequence analysis of 17 isolates from human and environmental sources revealed 13 clusters of strains and identified them as Burkholderia glathei Glade (BGC) bacteria. The taxonomic status of these clusters was examined by whole-genome sequence analysis, determination of the G+C content, whole-cell fatty acid analysis and biochemical characterization. The whole-genome sequence-based phylogeny was assessed using the Genome Blast Distance Phylogeny (GBDP) method and an extended multilocus sequence analysis (MLSA) approach. The results demonstrated that these 17 BGC isolates represented 13 novel Burkholderia species that could be distinguished by both genotypic and phenotypic characteristics. BGC strains exhibited a broad metabolic versatility and developed beneficial, symbiotic, and pathogenic interactions with different hosts. Our data also confirmed that there is no phylogenetic subdivision in the genus Burkholderia that distinguishes beneficial from pathogenic strains. We therefore propose to formally classify the 13 novel BGC Burkholderia species as Burkholderia arvi sp. nov. (type strain LMG 29317(T) = CCUG 68412(T)), Burkholderia hypogeia sp. nov. (type strain LMG 29322(T) = CCUG 68407(T)), Burkholderia ptereochthonis sp. nov. (type strain LMG 29326(T) = CCUG 68403(T)), Burkholderia glebae sp. nov. (type strain LMG 29325(T) = CCUG 68404(T)), Burkholderia pedi sp. nov. (type strain LMG 29323(T) = CCUG 68406(T)), Burkholderia arationis sp. nov. (type strain LMG 29324(T) = CCUG 68405(T)), Burkholderia fortuita sp. nov. (type strain LMG 29320(T) = CCUG 68409(T)), Burkholderia temeraria sp. nov. (type strain LMG 29319(T) = CCUG 68410(T)), Burkholderia calidae sp. nov. (type strain LMG 29321(T) = CCUG 68408(T)), Burkholderia concitans sp. nov. (type strain LMG 29315(T) = CCUG 68414(T)), Burkholderia turbans sp. nov. (type strain LMG 29316(T) = CCUG 68413(T)), Burkholderia catudaia sp. nov. (type strain LMG 29318(T) = CCUG 68411(T)) and Burkholderia peredens sp. nov. (type strain LMG 29314(T) = CCUG 68415(T)). Furthermore, we present emended descriptions of the species Burkholderia sordidicola, Burkholderia zhejlangensis and Burkholderia grimmiae. The GenBank/EMBUDDBJ accession numbers for the 16S rRNA and gyrB gene sequences determined in this study are LT158612-LT158624 and LT158625-LT1158641, respectively

Use of Whole Genome Phylogeny and Comparisons in the Development of a Multiplex-PCR Assay to Identify Sequence Type 36 Vibrio parahaemolyticus

Vibrio parahaemolyticus sequence type (ST) 36 strains that are native to the Pacific Ocean have recently caused multi-state outbreaks of gastroenteritis linked to shellfish harvested from the Atlantic Ocean. Whole genome comparisons of 295 genomes of V. parahaemolyticus, including several traced to northeastern US sources, were used to identify diagnostic loci: one putatively encoding an endonuclease (prp), and two others potentially conferring O-antigenic properties (cps and flp). The combination of all three loci was present only in one clade of closely-related strains, of ST36, ST59 and one additional unknown sequence type. However, each locus was also identified outside this clade, with prp and flp occurring in only two non-clade isolates, and cps in four. Based on the distribution of these loci in sequenced genomes, prp could identify clade strains with \u3e99% accuracy, but the addition of one more locus would increase accuracy to 100%. Oligonucleotide primers targeting prp and cps were combined in a multiplex PCR method that defines species using the tlh locus, and determines presence of both the tdh and trh hemolysin-encoding genes which are also present in ST36. Application of the method in vitro to a collection of 94 clinical isolates collected over a four year period in three Northeastern US, and 87 environmental isolates, revealed the prp and cps amplicons were only detected in clinical isolates identified as belonging to the ST36-clade, and in no environmental isolates from the region. The assay should improve detection and surveillance, thereby reducing infections

Susceptibility of Caenorhabditis elegans to Burkholderia Infection Depends on Prior Diet and Secreted Bacterial Attractants

The nematode Caenorhabditis elegans may be killed by certain pathogenic bacteria and thus is a model organism for studying interactions between bacteria and animal hosts. However, growing nematodes on prey bacteria may influence their susceptibility to potential pathogens. A method of axenic nematode culture was developed to isolate and quantify interactions between C. elegans and potentially pathogenic strains of the Burkholderia cepacia complex. Studying these dynamics in liquid solution rather than on agar surfaces minimized nematode avoidance behavior and resolved more differences among isolates. Most isolates of B. cenocepacia, B. ambifaria and B. cepacia caused 60–80% mortality of nematodes after 7 days, whereas isolates of B. multivorans caused less mortality (<25%) and supported nematode reproduction. However, some B. cenocepacia isolates recovered from chronic infections were much less virulent (5–28% mortality). As predicted, prior diet altered the outcome of interactions between nematodes and bacteria. When given the choice between Burkholderia and E. coli as prey on agar, axenically raised nematodes initially preferred most lethal Burkholderia isolates to E. coli as a food source, but this was not the case for nematodes fed E. coli, which avoided toxic Burkholderia. This food preference was associated with the cell-free supernatant and thus secreted compounds likely mediated bacterial-nematode interactions. This model, which isolates interactions between bacteria and nematodes from the effects of prior feeding, demonstrates that bacteria can influence nematode behavior and their susceptibility to pathogens

Why Genes Evolve Faster on Secondary Chromosomes in Bacteria

In bacterial genomes composed of more than one chromosome, one replicon is typically larger, harbors more essential genes than the others, and is considered primary. The greater variability of secondary chromosomes among related taxa has led to the theory that they serve as an accessory genome for specific niches or conditions. By this rationale, purifying selection should be weaker on genes on secondary chromosomes because of their reduced necessity or usage. To test this hypothesis we selected bacterial genomes composed of multiple chromosomes from two genera, Burkholderia and Vibrio, and quantified the evolutionary rates (dN and dS) of all orthologs within each genus. Both evolutionary rate parameters were faster among orthologs found on secondary chromosomes than those on the primary chromosome. Further, in every bacterial genome with multiple chromosomes that we studied, genes on secondary chromosomes exhibited significantly weaker codon usage bias than those on primary chromosomes. Faster evolution and reduced codon bias could in turn result from global effects of chromosome position, as genes on secondary chromosomes experience reduced dosage and expression due to their delayed replication, or selection on specific gene attributes. These alternatives were evaluated using orthologs common to genomes with multiple chromosomes and genomes with single chromosomes. Analysis of these ortholog sets suggested that inherently fast-evolving genes tend to be sorted to secondary chromosomes when they arise; however, prolonged evolution on a secondary chromosome further accelerated substitution rates. In summary, secondary chromosomes in bacteria are evolutionary test beds where genes are weakly preserved and evolve more rapidly, likely because they are used less frequently